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MULTIOBJECTIVE FRACTIONAL VARIATIONAL CALCULUS 
IN TERMS OF A COMBINED CAPUTO DERIVATIVE* 
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Abstract. The study of fractional variational problems in terms of a combined fractional Caputo 
derivative is introduced. Necessary optimality conditions of Euler— Lagrange type for the basic, 
isoperimetric, and Lagrange variational problems are proved, as well as transversality and sufficient 
optimality conditions. This allows to obtain necessary and sufficient Pareto optimality conditions 
for multiobjective fractional variational problems. 
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1. Introduction. There is an increasing interest in the study of dynamic sys- 
tems of fractional (where "fractional" actually means "non- integer" ) order. Extending 
derivatives and integrals from integer to non-integer order has a firm and longstand- 
ing theoretical foundation. Leibniz mentioned this concept in a letter to L'Hopital 
over three hundred years ago. Following L'Hopital's and Leibniz's first inquisition, 
fractional calculus was primarily a study reserved to the best minds in mathematics. 
Fourier, Euler, and Laplace are among the many that contributed to the development 
of fractional calculus. Along the history, many found, using their own notation and 
methodology, definitions that fit the concept of a non-integer order integral or deriva- 
tive. The most famous of these definitions among mathematicians, that have been 
popularized in the literature of fractional calculus, are the ones of Riemann-Liouville 
and Grunwald-Letnikov. On the other hand, the most intriguing and useful applica- 
tions of fractional derivatives and integrals in engineering and science have been found 
in the past one hundred years. In some cases, the mathematical notions evolved in or- 
der to better meet the requirements of physical reality. The best example of this is the 
Caputo fractional derivative, nowadays the most popular fractional operator among 
engineers and applied scientists, obtained by reformulating the "classical" definition 
of Riemann-Liouville fractional derivative in order to be possible to solve fractional 
initial value problems with standard initial conditions |33]. Particularly in the last 
decade of the XX century, numerous applications and physical manifestations of frac- 
tional calculus have been found. Fractional differentiation is nowadays recognized 
as a good tool in various different fields: physics, signal processing, fluid mechanics, 
viscoelasticity, mathematical biology, electrochemistry, chemistry, economics, engi- 
neering, and control theory (see, e.g., [9l[l0l[Ta[l9l[2Ql[23l[25l|30l[37l|4l|4^ 

The fractional calculus of variations was born in 1996 with the work of Riewe 
[321110], and is nowadays a subject under strong current research (see PlISllTllSl fTTlfT^ 
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[T5HT71I55] and references therein). The fractional calculus of variations extends the 
classical variational calculus by considering fractional derivatives into the variational 
integrals to be extremized. This occurs naturally in many problems of physics and 
mechanics, in order to provide more accurate models of physical phenomena (see, 
e.g., lUI^). The aims of this paper are twofold. Firstly, we extend the notion of 
Caputo fractional derivative to the fractional derivative "^I?"' , which is a convex 
combination of the left Caputo fractional derivative of order a and the right Caputo 
fractional derivative of order /?. This idea goes back at least as far as [22,, where based 
on the Riemann-Liouville fractional derivatives, the symmetric fractional derivative 
was introduced. Klimek's approach [22 is obtained in our framework as a particular 
case, by choosing parameter 7 to be 1/2. Although the symmetric fractional deriva- 
tive of Riemann-Liouville introduced by Klimek is a useful tool in the description 
of some nonconservative models, this type of differentiation does not seems suitable 
for all kinds of variational problems. Indeed, the hypothesis that admissible trajecto- 
ries y have continuous symmetric fractional derivatives implies that y{a) = y{b) = 
(cf. [41]). Therefore, the advantage of the fractional Caputo-type derivative '-"D'^'^ 
here introduced lies in the fact that using this derivative we can describe a more 
general class of variational problems. It is also worth pointing out that the fractional 
derivative '-'D"'^ allows to generalize the results presented in [3]. Our second aim 
is to introduce the subject of multiobjective fractional variational problems. This 
seems to be a completely open area of research, never considered before in the lit- 
erature. Knowing the importance and relevance of multiobjective problems of the 
calculus of variations in physics, engineering, and economics (see [13,26,27,33,43,46] 
and the references given there), and the usefulness of fractional variational problems, 
we trust that the results now obtained will open interesting possibilities for future 
research. Main results of the the paper provide methods for identifying Pareto op- 
timal solutions. Necessary and sufficient Pareto optimality conditions are obtained 
by converting a multiobjective fractional variational problem into a single or a family 
of single fractional variational problems with an auxiliary scalar functional, possibly 
depending on some parameters. 

The paper is organized as follows. Section [2] presents some preliminaries on frac- 
tional calculus, essentially to fix notations. In Section [3] we introduce the fractional 
derivative D"'^ and provide the necessary concepts and results needed in the sequel. 
Our main results are stated and proved in Section |4] and Section [5] The fractional 
variational problems under our consideration are formulated in terms of the fractional 
derivative ^D"^'^ . We discuss the fundamental concepts of a variational calculus such 
as the the Euler-Lagrange equations for the elementary fSubsection l4.1|) . isoperimetric 
(Subsection 14. 3p . and Lagrange (Subsection 14. 4|) problems, as well as sufficient opti- 
mality fSubsection l4.5p and transversality (Subsection l4.2|) conditions. Section[5]deals 
with the multiobjective fractional variational calculus. We present Pareto optimality 
conditions (Subsection 15. ip and examples illustrating our results (Subsection 15.21) . 

2. Fractional calculus. In this section we review the necessary definitions and 
facts from fractional calculus. For more on the subject we refer the reader to the 
books pni^l5glH^ . Let / G Li([a, &]) and < a < 1. We begin by defining the left 
and the right Riemann-Liouville Fractional Integrals (RLFI) of order a of a function 
/. The left RLFI is given by 




X e [a, b], 



(2.1) 
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and the right RLFI by 



1 '•^ 



^Ibfi^) -Ffr / (i " ^r~'fm, X G [a, b], (2.2) 
where r(-) represents the Gamma function, i.e., 

r(z) / e-^e-'dt. 



Moreover, Jxf = xl^f = / if / is a continuous function. The left and the right 
Riemann-Liouville derivatives are defined with the help of the respective fractional 
integrals. The left Riemann-Liouville Fractional Derivative (RLFD) is given by 

aD'^Ji^) ■■= Jj"^ - = j-jl-'-f{x), X e [a, 6], (2.3) 

and the right RLFD by 

-1 d /"V . , ( d 



Let / € AC{[a,b]), where AC{[a,b]) represents the space of absolutely continuous 
functions on [a,b]. Then the Caputo fractional derivatives are defined as follows: the 
left Caputo Fractional Derivative (CFD) by 

^D^fix) := \ Jjx ^ t)-"^^f{t)dt ^ Jl--±f{x), a;e[a,6], (2.5) 

and the right CFD by 

'.Dtm ^ £(t-.)-"|/(.).t . Jl-- /(.), . e (2.6) 



where a is the order of the derivative. The operators (|2.ip - (|2.6p are obviously linear. 
We now present the rule of fractional integration by parts for RLFI (see, e.g., [M])- 
Let < a < 1, p > 1, > 1, and l/p+l/q < 1 + a. If g e Lp{[a, b]) and / e Lq{[a, 6]), 
then 

g{x)al^f{x)dx = / f{x)xl^g{x)dx. (2.7) 

J a 

In the discussion to follow, we will also need the following formulae for fractional 
integrations by parts: 



g{x)^^D'if{x)dx^ f{xW,-'^g{x)l^^+ / f{x)xD'ig{x)dx, 

6 (2.8) 
g{x)<:j^D^!{x)dx = -f{x)Jl~^g{x)^;:^^ + f f{x)aD'^^g{x)dx. 



They can be derived using equations (I2.3l) - (|2.6p . the identity (|2.7p . and performing 
integration by parts. 
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3. The fractional operator '^D"'^ . Let a, /3 G (0, 1) and 7 e [0, 1] . We define 



the fractional derivative operator ^D"'^ by 

^D-/:=^'^D: + {1~j)^'D^,, (3.1) 
wfiicfi acts on / S AC {[a, b]) in the expected way: 

^D;-''fix) = ^'^D^fix) + (1 - 7)?^f /(x). 

Note that ^D'^'^f{x) = ^D^f{x) and ^D'^'^f{x) = ^D^f{x). The operator (IXTj) 
is obviously linear. Using equations (12.81) it is easy to derive the following rule of 
fractional integration by parts for D"'^: 

' g{x) ^D;'^f{x)dx = 7 [f{xUl-^g{x)\ 

+ (1 - 7) [-f{x)Jl-^g{x)] + f f{x)Dt"^gix)dx, (3.2) 

J a 

where i^fl"^ := (1 - j)aD^^ + -f^D^. Let e N and f = [A, . . . , /jv] : [a, b] ^ 
with /, e AC{[a,b]), i = 1, . . . , TV; a, /3, 7 G with a,,/3, G (0,1) and 7, G [0,1], 
i = 1,...,7V. Then, 

^i^^'V) pi^^;'^^/i(:r),...,^i?:;~'^"/iv(x)] . 

Let D denote the set of all functions y : [a, b] — such that '-^ D'^'^y exists and is 
continuous on the interval [a, b\. We endow D with the following norm: 

||y||i,oo := max ||y(x)|| + max ||'^i:'"''^y(x)||, 

a<x<b a<x<b 

where \\ ■ || is a norm in M^. Along the work we denote by diK, i = 1, . . . , M [M G N), 
the partial derivative of function K : M*'^ — )■ R with respect to its ith argument. Let 
A G R''. For simplicity of notation we introduce the operators [y]"''^ and A{y}"''^ by 

[y]^-^ix) := [xM^),'''D"/yi^)) , 
x{yr/{x) := [x,y{x),''D;'\{x),Xu...,Xr) ■ 

4. Calculus of variations via ^B"/. We are concerned with the problem 
of finding the minimum of a functional j' : V ^ M., where I? is a subset of D. 
The formulation of a problem of the calculus of variations requires two steps: the 
specification of a performance criterion, and the statement of physical constraints 
that should be satisfied. The performance criterion J7, also called cost functional 
(or objective), must be specified for evaluating quantitatively the performance of the 
system under study. We consider the following cost: 

Jiy)= [\[y]^^^ix)dx, 
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where x G [a, h] is the independent variable; y{x) € is a real vector variable, the 
functions y are generally called trajectories or curves; ^D"'^y{x) G stands for 
the fractional derivative of y{x); and L £ C^([a, h] x M^^; M) is the Lagrangian. 

Enforcing constraints in the optimization problem reduces the set of candidate 
functions and leads to the following definition. 

Definition 4.1. A trajectory y G D is said to he an admissible trajectory, 
provided it satisfies all the constraints of the problem along the interval [a, h]. The set 
of admissible trajectories is defined as T) :— {y ETi : y is admissible} . 

We now define what is meant by a minimum of J7 on P. 

Definition 4.2. A trajectory y G I? is said to be a local minimizer for J on T>, 
if there exists 6 > such that J{y) < J{y) for ally with \\y — y||i,oo < 5. 

The concept of variation of a functional is central to the solution of problems of 
the calculus of variations. 

Definition 4.3. The first variation of J at y ^ T) in the direction h G D is 
defined as 

SJ{y; h) := hm = -r^Jiy + eh) 

provided the limit exists. 

Definition 4.4. A direction h G D. h ^ 0, is said to be an admissible variation 
for J aty €V if 

(i) '^i7(y; h) exists; and 

(ii) y + eh G I? for all sufficiently small e. 

The following well known result offers a necessary optimality condition for the 
problems of the calculus of variations, based on the concept of variations. 

Theorem 4.5 (see, e.g.. Proposition 5.5 of [45]). Let J be a functional defined 
on T>. Suppose that y is a local minimizer for J on T>. Then, 5J{y; h) = for each 
admissible variation h. at y. 

4.1. Elementary problem of the '-'D""^ fractional calculus of variations. 

Let us begin with the following fundamental problem: 

Jiy)^ L[y]'^'^ix)dx^mm (4.1) 

J a 

over all y G D satisfying the boundary conditions 

y(a)=y^ y(fe) = y^ (4.2) 

where y°',y^ G are given. The next theorem gives the fractional Euler-Lagrange 
equation for the problem (|4.1[) " (|4.2p . 

Theorem 4.6. Let y = (yi, ...,?/ at) be a local minimizer to problem ()4.1|) - ()4.2|) . 

Then, y satisfies the following system of N fractional Euler-Lagrange equations: 

my]y'^i^) + i^fr;'_";-^5jv+.L[y]^'^(x) = O, z = 2,...N + 1, (4.3) 
for all X G [a, b] . 

Proof. Suppose that y is a local minimizer for J'. Let h be an arbitrary admissible 
variation for problem (|4.ip - (|4.2[) . i.e., hi{a) — hi{h) — 0, i — 1, . . . , N . Based on the 
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differentiability properties of L and Theorem 14. 5[ a necessary condition for y to be a 
local minimizer is given by 



d_ 

Ye 



J{y + eh) 



0, 



e=0 



that is 

l-b 



'N+l N+1 
. 1=2 i=2 



dx = 0. 



Using formulae 
function, we get 



(4.4) 

of integration by parts in the second term of the integrand 



■jv+i 



i=2 



hi^i{x)dx 



+ 7 



'N+l 



(1-7) 



'N+l 



x—b 



= 0. (4.5) 



Since hi{a) — hi(b) — 0, i — 1, . . . , N , hy the fundamental lemma of the calculus of 
variations we deduce that 



dMy]"^^{x) + Dt-^\^;-'^M+^L[y]^^^ix) = 0, z 



.,N + 1, 



Cn/3 

b 



for all X G [a,b]. U 

Observe that if a and /3 go to 1, then ^ D" can be replaced with ^ and ^ D[ 
with — ^ (see [3H]). Thus, if 7 = 1 or 7 = 0, then for a, 13 1 we obtain a 
corresponding result in the classical context of the calculus of variations (see, e.g., [45) 
Proposition 6.1]). 

4.2. Fractional transversality conditions. Let I G {1,. . . ,N}. Assume that 
y(a) = y°, yi{b) = y^, i — 1, . . . ,N , i ^ I, but yi{b) is free. Then, hi{b) is free and 
by equations (|4.3I) and (|4.5I) we obtain 



j.ll-'''dN+i+iL[y]^^^{x) -{l-j)j',-^'dN+i+iL[yr'0ix) 



0, 



x — b 



where a,/3,7 G K. Let us consider now the case when y(a) — y°, yi{b) = y^, 
i = 1, . . . , N , i I, and yi{b) is free but restricted by a terminal condition yi{b) < y\. 
Then, in the optimal solution y, we have two possible types of outcome: yi{b) < y\ or 
yi{b) = yf. If yi{b) < yf, then there are admissible neighboring paths with terminal 
value both above and below yi{b), so that hi{b) can take either sign. Therefore, the 
transversality conditions is 



-f.jl-^'dN+i+iLiyr^f'ix) -{l--f)al',-^'dN+i+iL[yr-^{x) 



= 



x — b 



(4.6) 



for yi{b) < y^. The other outcome yi{b) = y^ only admits the neighboring paths with 
terminal value yi{b) < yi{b). Assuming, without loss of generality, that hi{b) > 0, this 
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means that e < 0. Hence, the transversahty condition, which has it root in the first 
order condition (|4.4p . must be changed to an inequahty. For a minimization problem, 
the < type of inequahty is cahed for, and we obtain 



7./^"'aw+i+ii[y]"^^(:E) -(1 -7)a/^^'9w+i+iL[y]"'^(x) 



< 



(4.7) 



for iji{b) = yf. Combining (|4.6p and (14. 7p . we may write the following transversahty 
condition for a minimization problem: 



x—h 



(yiib) - y'l) 7./r"'5jv+/+ii[y]"'^(a;) -(1 - ^)Jl-('^dN+i+iL[yr^-^{x) 



0. 



4.3. The *^D"'^ fractional isoperimetric problem. Let us consider now the 
isoperimetric problem that consists of minimizing (14. ip over all y G D satisfying r 
isoperimetric constraints 



(4.8) 



where £ C^([a, h] x R^^; M), j = 1, . . . , r, and boundary conditions 

Necessary optimality conditions for isoperimetric problems can be obtained by 
the following general theorem. 

Theorem 4.7 (see, e.g.. Theorem 2 of [TH] on p. 91). LetJ.G^,..., g'' be func- 
tionals defined in a neighborhood of y and having continuous first variations in this 
neighborhood. Suppose that y is a local minimizer of (|4.ip subject to the boundary 
conditions (14. 2p and the isoperimetric constrains (j4.8p . Assume that there are func- 
tions h^, . . . , h'' G D such that the matrix A — (0^1), a^i := 5Q^{y\ h'), has maximal 
rank r. Then there exist constants Ai, . . . , G K such that the functional 



satisfies 



5^(y;h) =0 



(4.9) 



for all heD 

Suppose now that assumptions of Theorem 14.71 hold. Then, equation (|4.9p is 
fulfilled for every h G D. Let us consider function h such that h(a) = h(6) = 0. 

Then, we have 

d 



= <5^(y;h) = ^(y + eh)|,=o 
oe 

'N+l N+1 

E d^,F,{y}^^P{x)h,^,{x) + E ^N+^Fx{y}^^^{xfD^::l^^^-^h,^,{x) 

. i=2 1=2 



dx. 



where the function F : [a, b] x 



p27V 



is defined by 



F^{yr^'\x) L[yr^^l\x) ^.G-' [y]?^^(^ 
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On account of the above, and similarly in spirit to the proof of Theorem l4.6[ we obtain 
d,xF{y}^-^{x) + Dt-,]•T"9N+^xF{y}^^^ix) =0, * = 2, . . .iV + 1. (4.10) 

Therefore, we have the following necessary optimality condition for the fractional 
isoperimetric problems: 

Theorem 4.8. Let assumptions of Theorem \4.7 hold. If y is a local minimizer 
to the isoperimetric problem given by (j4.ip . (j4.2p and (14.81) . then y satisfies the system 
of N fractional Euler- Lagrange equations (|4.10p for all x G [a, b] . 

Suppose now that constraints (|4.8p are characterized by inequalities 



Q'iy) = / G^[yp^ix)dx < Z„ J = 1, . . . ,r. 

J a 

In this case we can set 

(^Gny]^^^{x) ~ ^) dx + j\<P,{x)fdx = 0, 

j = 1, . . . , r, where (pj have the some continuity properties as yi. Therefore, we obtain 
the following problem: 

pb 

J{y)'=^ I Hx,y{x),'^ d"' y{x),(t)i{x),...,<j)r{x))dx — ^ min 



dx^O, j = l,...,r. 



subject to r isoperimetric constraints 



and boundary conditions (|4.2p . Assuming that assumptions of Theorem 14.81 are satis- 
fied, we conclude that there exist constants Xj € M., j — 1, . . . ,r, for which the system 
of equations 

£'fi"^V"r'^w+''^(^' yix)i'^ F)"'^y{x), Ai, . . . , Ar, 0i(a;), . . . , (t)rix)) 

+ d,F{x, yix), ^i?"'^y(x), Ai, . . . , A,, M^), ■ • ■ , M^)) = 0, (4.11) 
i = 2,...,7V + l, F = L + E;=i a, (G^' - ^ + ^^,) and 

A,0,(x)=O, j = l,...,r, (4.12) 

hold for all x G [a,b]. Note that it is enough to assume that the regularity condition 
holds for the constraints which are active at the local minimizer y (constraint Q'' is 
active at y if Q'^{y) = h)- Indeed, suppose that I < r constrains, say Q^, . . . ,Q'' for 
simplicity, are active at the local minimizer y, and there are functions h^, . . . , h' £ D 
such that the matrix 

B^ibkj), :=(5g^-(y;h-'), k,j^l,...,l<r 

has maximal rank I. Since the inequality constraints . .. are inactive, the 

condition (|4.12l) is trivially satisfied by taking Aj+i = • • ■ = A^ = 0. On the other 
hand, since the inequality constraints G^,---,G^ are active and satisfy a regularity 
condition at y, the conclusion that there exist constants Xj G K, j — 1, . . . , r, such 
that (|4.1ip holds follow from Theorem 14.81 Moreover, (|4.12l) is trivially satisfied for 
j = l,...,L 
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4.4. The *^D"'^ fractional Lagrange problem. Let us consider tlie following 
Lagrange problem, which consists of minimizing (j4.ip over all y G D satisfying r 



independent constraints (r < N) 



Gny]"^^{x) = 0, .7 = 1,.. 



(4.13) 

and boundary conditions (|4.2p . In mechanics, constraints of type (|4.13p are called 
nonholonomic. By the independence of the r constraints G C^{[a,b] x R^^;M) it 
is meant that it should exist a nonvanishing Jacobian determinant of order r, such as 

■a? — — — — - — !■ 7^0. Of course, any r of p, , 7 = + 2, 2A^ + 1, can be used, not 
necessarily the first r. 

Theorem 4.9. A junction y which is a solution to problem (|4.H) . (|4.2|) subject to 
r independent constraints (r < N ) ()4.13p satisfies, for suitably chosen functions Xj, 
j = 1, . . . , r, the system of N fractional Euler-Lagrange equations 

^,F[y,\]^^^ix)+Dt-^\'";-'^N+^F[y,X]^^^{x)^0, x e [a,b], z^2,...,N + l, 

where F[y , X]'^'^ (x) = L[yp0 (x) + A,(2:)GJ'[y]^'^(x). 

Proof Suppose that y — (yi, . . . , 2/Ar) is the solution to problem defined by 
()4.ip . ()4.2p . and ()4.13p . Let h = {hi, . . . , Hn) be an arbitrary admissible variation, 
i.e., /i,(a) = = 0, i = 1, . . . , A^, and G^[y + shp^{x) = 0, j = 1, . . . ,r, where 

e e M is a small parameter. Because y = {yi, . . . , yjsr) is a solution to problem defined 
by (|4IT]),(I43), and (|4l3l) . it follows that 



|j(y + .h) 



0, 



e=0 



that is. 



'N+l 



N+1 



^ dMy]^-^ix)h,.i{x) + J2 dN+rL[y]^-^ixfD^::^^^^-^h,.iix) 



and, for j = I, . . . ,r. 



dx — 0, 
(4.14) 



N+1 N+1 

d,G'[y]^'^{x)h,_i{x) + J2 ^N+^G^[y]^'^{xfD^^:l'^^-'h,_l{x) = 0. (4.15) 

Multiplying the jth equation of the system (|4.15p by the unspecified function Xj{x), 
for all j = 1, . . . , r, integrating with respect to x, and adding the left-hand sides (all 
equal to zero for any choice of the Xj ) to the integrand of (|4.14p , we obtain 

N+1 



J2 my]^^^{x)+J2>^Ji^)^^Gny]^'^{x) K^i^x) 

i=2 \ i=l ) 

V+1 I ^ \ 

+ f a^+.i[y]^;'''(x) + ^ A,(x)ajv+.G-'[y]^^'^(x) j (^i?:;;_--ft-/j.„i(x)) 

N+1 N+1 

J2 d,F[y,Xp^{x)h,.i{x) + J2 ^N+^F[y,Xp^{x){^D^::^'^'-^h.,.lix)) 



i=2 



dx 



dx 



0, 
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where F[y,Xp^{x) = L[yp^{x) + ^-^^ X{x)&[y]^'^ {x). Integrating by parts, 



'N+l 



J2 ^^F[y, \]^^P{x) + Dt^:^l-'^N+^F[y, XT^^^ {x) 



i=2 



h.,-i{x)dx 0. (4.16) 



Because of (|4.15p , we cannot regard the A'' functions hi, . . . ,hN as being free for arbi- 
trary choice. There is a subset of r of these functions whose assignment is restricted 
by the assignment of the remaining (N — r). We can assume, without loss of gener- 
ality, that hi, . . . ,hr are the functions of the set whose dependence upon the choice 
of the arbitrary hr+i, . . . ,hN is governed by (|4.15|) . We now assign the functions 
Xi, . . . , Xr to be the set of r functions that make vanish (for all x between a and b) 
the coefficients oi hi, . . . , hr in the integrand of (|4.16l) . That is, Ai, . . . , are chosen 
so as to satisfy 

d,F[y, X]^-^ix) + Dt-^;^l-'^N+^F[y, X^^-'^ix) =0, * = 2, . . . ,r + 1, x<-[a,h]. 



With this choice (|4.16p gives 



(4.17) 



J a 



N+l 



^^F[y, >^r^'^{^) + Dt^;^-'^N+^F[y, XT^^^ix) 



.i=r+2 



hi-i{x)dx = 0. 



Since the functions /ir+i, ■ • ■ , ^tv are arbitrary, we may employ the fundamental lemma 
of the calculus of variations to conclude that 



a,F[y, Xr/ix) + DP-^;^:-'dr,+,F[y,Xr^'^{x) = 



(4.18) 



i = r + 1, . . . , + 1, for all X G [a, b]. □ 

Remark 4.10. In order to determine the {N -f r) unknown functions yi, . . . , y„, 
Ai, . . . , Xr, we must consider the system of (N + r) equations, consisting of (j4.13p . 
()4.17p . and (|4.18p . together with the 2N boundary conditions ()4.2p . 

Assume now that the constraints, instead of (|4.13p . are characterized by inequal- 
ities: 

G^y]^'^ix)<0, j = l,...,r. 

In this case we can set 

G^[y];;^^(a:) + (0,(.T))2=O, J = l,...,r, 

where have the some continuity properties as yi. Therefore, we obtain the following 
problem: 



Jiy) ^ / Lix,y{x),^D"'^y{x),(t)i{x),...,(j)r{x))dx - 

J a 

subject to r independent constraints (r < N) 

G^[y];;^^(a:) + (0,(:r))2=O, J^l,...,r, 



(4.19) 



(4.20) 



and boundary conditions (|4.2p . Applying Theorem 14.91 we get the following result. 

Theorem 4.11. A set of functions yi, . . . ,yN , (t^i, ■ ■ ■ , (t^r, which is a solution 
to problem (j4.19p - (l4.20p . satisfies, for suitably chosen Xj , j — 1, . . . ,r, the following 
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system of equations: 

D^l-^'^'^^-'dN+iF^x, y{x),'^'D""'^y{x), Ai(x), . . . , Xr{x), (j)i{x), (I)r{x)) 

+ dtF{x, y{x), ^L'"'^y(a;), Xi{x), \r{x), . . . , = 0, 

i = 2,...,N+l, where F = i + Ej=i (G^ + ^i). and Xj{x)(j)j{x) = 0, j = 1, . . . ,r, 
hold for all x G [a,b]. 

4.5. Sufficient condition of optimality. In this section we provide sufficient 
optimality conditions for the elementary and the isoperimetric problem of the D"'^ 
fractional calculus of variations. Similarly to what happens in the classical calculus 
of variations, some conditions of convexity are in order. 

Definition 4.12. Given a function f e C^([a, b] xR^^; M), we say that f{x, y, v) 
is jointly convex in (y, v), if 

N+l N+1 

/(x,y + y°,v + v") -/(a;,y,v) > ^ d,f{x,y,^)yl, + ^ dN+^f{x,y,^)vl, 

for all (a;,y,v),(a;,y +y",v + v") E [a,b] x R"^^ . 

Theorem 4.13. Let L{x,y,v) be jointly convex in (y,v). If y satisfies the 
system of N fractional Euler-Lagrange equations (|4.3p , then y is a global minimizer 
to problem 

Proof. The proof is similar to the proof of Theorem 3.3 in [2^. □ 
Theorem 4.14. Let F{x,y,v,X) = L{x,y,v) - Y^^j^i^jG^ {x,y,w) be jointly 
convex in (y,v), for some constants Xj G M., j — 1, . . . ,r. If y^ satisfies the sys- 
tem of N fractional Euler-Lagrange equations (|4.10l) . then y^ is a minimizer to the 
isoperimetric problem defined by (j4.ip . (j4.2p and (j4.8p . 

Proof. By Theorem 14.131 y° minimizes F^^y'^'^'^ [x) dx. That is, for all func- 
tions satisfying condition (|4.2p we have 

/' L[yr^'\x) dx-J2 h f G=[yr^'^{x) dx 

> f L[yX'^{x)dx-j2h f G^[yX^'\x)dx. 

J a J — 1 



Restricting to the isoperimetric constraints (j4.8p . we obtain that 

/ L[y]^^^(a;) rfx - ^ X^lj dx > L[y°]'^^'^ (x) dx '^jh dx. 

Ja Ja j^j^ 

Therefore, 

L[y]^-'\x)dx> f L[yX'^{x)dx 



as desired. □ 

Choosing r = 1 in Theorem 14. 141 one can easily obtain ^ Theorem 3.10] 
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5. Multiobjective fractional optimization. Multiobjective optimization is a 
natural extension of the traditional optimization of a single-objective function. If the 
objective functions are commensurate, minimizing one-objective function minimizes 
all criteria and the problem can be solved using tradicional optimization techniques. 
However, if the objective functions are incommensurate, or competing, then the min- 
imization of one objective function requires a compromise in another objective. Here 
we consider multiobjective fractional variational problems with a finite number d > 1 
of objective (cost) functionals 

(J^(y),...,j'^(y)) = l^j\\yr^^^{x)dx,..., j\%r^^^{x)d:^ ^min (5.1) 

subject to the boundary conditions 

y(a) = y^ y{b) = Y\ (5.2) 

y°' £ M^, and r (r < N) independent constraints 

G^[yr/\x)<0, j = l,...,r, (5.3) 

where L\ G C^i[a, b] x R^w. r)^ i = 1, . . . , d, j = 1, . . . , r. We would like to find 
a function y G D, satisfying constraints (15. 2p and (j5.3p . that renders the minimum 
value to each functional J^, i = 1, . . . ,d, simultaneously. The competition between 
objectives gives rise to the necessity of distinguish between the difference of multiob- 
jective optimization and traditional single-objective optimization. Competition causes 
the lack of complete order for multiobjective optimization problems. The concept of 
Pareto optimality is therefore used to characterize a solution to the multiobjective 
optimization problem. For the usefulness of variational analysis and Pareto optimal 
allocations in welfare economics, we refer the reader to [32]. We define 

£ := {y G D : y satisyies conditions (|5.2p and (|5.3p }. 

Definition 5.1. A function y € £ is called a Pareto optimal solution to problem 
(|5.ip - (|5.3p if does not exist y € £ with 



Vie {!,..., d}:r{y)<r{y) A 3z G {1, . . . , d} : J^(y) < J'(y). 



Definition 15.11 introduces the notion of global Pareto optimality. Another impor- 
tant concept is the one of local Pareto optimality. 

Definition 5.2. A function y ^ £ is called a local Pareto optimal solution 
to problem (|5.ip - (l5.3p if there exists 5 > Q for which does not exist y d £ with 
\\y ~ y||i,oo < S and 

yie{l,...,d}:r{y)<r{y) A 3te{l,...,d}:r{y)<T{y). 



Naturally, any global Pareto optimal solution is locally Pareto optimal. For 
enhanced notions of Pareto optimality of constrained multiobjective problems, the 
reader is referred to [6 . 
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5.1. Fractional Pareto optimality conditions. We obtain a sufficient condi- 
tion for Pareto optimality by modifying the original multiobjective fractional problem 
(|5.ip - (|5.3p into the following weighting problem: 



d 

L'[y]':;-'^{x)dx — ^min (5.4) 



subject to y G £, where Wi >0 for alH = 1, . . . , d, and X]i=i = ^■ 

Theorem 5.3. The solution of the weighting problem (|5.4p is Pareto optimal if 
the weighting coefficients are positive, that is, Wi > for all i — 1, . . . ,d. Moreover, 
the unique solution of the weighting problem (|5.4I) is Pareto optimal. 

Proof. Let y e f be a solution to problem (|5.4p with Wi > for alH — 1, . . . , d. 
Suppose that y is not Pareto optimal. Then, there exists y such that i7*(y) < J^iS) 
for all i — and J'-'(y) < J-'{y) for at least one j. Since Wi > for all 

i — 1, . . . ,d,we have X^iLi WiJ^y) < X^iLi WiJ^{y). This contradicts the minimality 
of y. Now, let y be the unique solution to (|5.4I) . If y is not Pareto optimal, then 
X^iLi WiJ^{y) < J2i=i ^!'^*(y)- This contradicts the uniqueness of y. □ 

Therefore, by varying the weights over the unit simplex {w = {wi, . . . , Wd) ■ Wi > 
= 1} ones obtains, in principle, different Pareto optimal solutions. The 
next theorem provides a necessary and sufficient condition for Pareto optimality. The 
result is analogous to that valid for the finite dimensional case (see, e.g.. Chapter 3.1 
and Chapter 3.3 of [31J. 

Theorem 5.4. A function y £ £ is Pareto optimal to problem (j5.ip - (|5.3p if and 
only if it is a solution to the scalar fractional variational problem 

V[y]'^'^{x)dx ^niin 

J a 

subject to y G 8 and 

''\^[y]"^^{x)dx< f V[y]^'P{x)dx, j = l,...,d, j ^ i, 



for each i — 1, . . . , d. 

Proof. Suppose that y is Pareto optimal. Then y £ Ck — {y E £ ■ JHy) < 
^"'(y)iJ — 1, • ■ • ,d,j 7^ fc} for all fc, so Ck ^0. If y does not minimize J^{y) on 
the constrained set Ck for some fc, then there exists y £ £ such that J^{y) < J^{y) 
and J-'{y) < J^'iy) for all j ^ fc. This contradicts the Pareto optimality of y. 
Now, suppose that y minimize each J7'^(y) on the constrained set Ck- If y is not 
Pareto optimal, then there exists y such that J^{y) < J^{y) for alH = 1, . . . , cZ and 
'^"' (y) < '^"' (y) for at least one j. This contradicts the minimality of y for J-'{y) on 
Cr □ 

Remark 5.5. For a function y £ £ to be Pareto optimal to problem (|5.ip ^ (|5.3p . 
it is necessary to be a solution to the fractional isoperimetric problems 



V[yr'^{x)dx 



subject to y £ £ and 



' V[y]"'\x)dx^ f U[y]'^'^{x)dx, j = l,...,d, 
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for all i — l,...,d. Therefore, necessary optimality conditions for the fractional 
isoperimetric problems (see Theorem \4-S^ ire also necessary for fractional Pareto 
optimality. 

5.2. Examples. We illustrate our results with two multiobjcctive fractional vari- 
ational problems. 

Example 1. Let y{x) — Ea{x°'), x G [0,1], where is the Mittag-Leffler 
function: 

oo ^. 

When a = 1, the Mittag-Leffler function is simply the exponential function: Ei{x) = 
e^. We note that the left Caputo fractional derivative of y is y (cf. f57|/, p. 98): 

^D^y{x)^yix). 

Consider the following multiobjcctive fractional variational problem (N = 1, 7 = 1, 
and d — 2): 

(jHy), J\y)) = (^j\oD^yi^)f dx, J^' y{x)^D:y{x) dx^ min (5.5) 
subject to 

y(0) = 0, yil)^ Ea.il). (5.6) 

Observe that y satisfies the necessary Pareto optimality conditions (see Remark\5_ 
Indeed, as shown in Example 1], y is a solution to the isoperimetric problem 

J\y)= f\s D^yi^))' dx ^ mm 
Jo 

subject to 

y{x)^D:y{x)dx^ / (yix)? dx. 
Jo 

Consider now the following fractional isoperimetric problem: 

J\y) = / y{x)'iD'^y{x) dx mm 
Jo 

subject to 

''\^DZy{x)f dx = l\^DZy{x)f dx. 



Let us apply Theorem |/^.7| The equality xDfy{x) — holds if and only if y{x) = 
d(l — x)"^^ with d G IR. (see \21\ Corollary 2.1]). Hence, y does not satisfy equation 
xDf{Q D'^y) = 0. The augmented function is 
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and the corresponding fractional Euler-Lagrange equation gives 

A solution to this equation is X — and y — y. Therefore, by Remark ] 5. 51 y = y is a 
candidate Pareto optimal solution to problem (|5.5l) - (|5.6p . 

Example 2. Consider the following multiobjective fractional variational problem: 

{J\y),j'{y)) = l^f^\{^,D'iy{x)^f{x)fdx,j'^\{^DZy{x)fdx^^u,in 

(5.7) 

subject to 

2/(0) ==0, 2/(1) -X, xeM, (5.8) 

where f is a fixed function. In this case we have N — 1, 7 = 1, and d — 2. By Theo- 
rem \5.3[ Pareto optimal solutions to problem (|5.7p - (|5.8p can be found by considering 
the family of problems 

wj^ ^{^D^y{x)-f{x))'dx + {l-w)J^ ^^D^y{x)dx^ mm (5.9) 
subject to 

2/(0) = 0, 2/(1) = X, xeR, (5.10) 

where w G [0,1]. Let us now fix w . By Theorem \4.6\ a solution to problem (|5.9p - (l5.10p 
satisfies the fractional Euler- Lagrange equation 

,D^{^D^y{x)~wf{x))=0. (5.11) 

Moreover, by Theorem \4-t3[ a solution to (j5.11l) is a global minimizer to problem 
(|5.9p - (|5.10p . Therefore, solving equation (|5.1ip for w € [0,1], we are able to obtain 
Pareto optimal solutions to problem (I5.7p - (|5.8p . In order to solve equation (j5.11l) . 
firstly we use Corollary 2.1 of f21f to get the following equation: 

^D^y{x)-wf{x)^d{l-x)°'-\ deR. (5.12) 

Equation (j5.12p needs to be solved numerically. We did numerical simulations using 
the MatLab solver fade for linear Fractional- Order Differential Equations (FODE) 
with constant coefficients, developed by Farshad Merrikh Bayat \29^ . The results for 
a = 1/2, f{x) — e^, and different values of the parameter w can be seen in Figure\57J\ 
Numerical results for different values of a show that when a — )■ 1 the fractional solu- 
tion converges to the solution of the classical problem of the calculus of variations. 

Acknowledgments. The authors would like to express their gratitude to Ivo 
Petras, for having called their attention to ^29^ as well as for helpful discussions on 
numerical aspects and available software packages for fractional differential equations. 
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